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ABSTRACT 



An investigation of stresses, arising from a thermal field, in 
electronic packages comprised of a ceramic chip attached to a 
composite board by a solder joint was undertaken. Two types of 
solder joint connections, a leaded device and an unleaded device, 
were analyzed. The stresses result from the mismatch of material 
properties of the tri-assembly system. Procedures are shown for the 
analyses of the two types of solder connections. The analyses 
showed that the thermoelastic stresses for the stiff unleaded 
connection are significantly greater, by a factor of 20 for a the 
particular application undertaken in this report, than the 

thermoelastic stresses in the flexible leaded connection. 



Introduction 



Electronic chip packages are comprised of several components with 
different material properties and geometries. One of the most basic 
electronic packages is the tri-assembly system consisting of a chip 
attached to a board by a solder connection. During the fabrication 
of this simple unit, the soldering process subjects the device to 
a thermal field. Since the three components of the device have 
different material properties, in particular their Young’s modulii, 
Poissons ratios, and thermal coefficients of expansion, the mismatch 
in these properties result in thermoelastic stresses. In addition to 
the thermal field arising from fabrication, thermal fields and hence 
thermoelastic stresses, result from the operation of the device itself. 

Not only is the mismatch in material properties important in the 
thermoelastic stresses that develop, but equally important is the 
manner in which the solder connection is made. For a given set 
of chip, board and solder properties, an important question ’what 
is the best way to make the solder connection?’. With regard to 
geometric configuration, there are two main types of connections, 
the leaded connection and the unleaded connection. In the leaded 
connection, the chip is attached to the board by a thin ’lead’ as 
shown in Figure 1. In this configuration, the chip is somewhat 
removed from the board its attached to. In contrast, in the 
unleaded connection, the chip is attached along its edges to the 
board substrate by a thin layer of solder. This unleaded joint 
arrangement is significantly ’stiffer’ than the leaded connection, and 
all other things being equal, one would expect that the stiffer 
unleaded joint would result in larger thermoelastic stresses. This 
investigation seeks to determine the relative magnitudes of the 
stresses from these two kinds of solder joints. It should be noted 
however, that the proximity of the chip to the board, in the 
unleaded joint, might result in a more uniform thermal field thereby 
reducing the thermal gradients and thus the thermal stresses. The 
question of the thermal field that results from each of the solder 
joints is not addressed in this investigation which only considers 
the effect of the same uniform thermal environment on each of 
these (leaded and unleaded) solder joints. 
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Analysis of leaded device 

Figure 1 shows the geometry and dimensions of the leaded 
device considered in this study which consists of a ceramic chip, 
lead wires and a circuit board. The cross-sectional dimension of the 
lead wire is given in Figure 2. 

Three general assumptions are made. First, the whole device is 
under a thermal cycling load of +65°c ~ -20°c and the lead wire has 
an initial strain from soldering process. The device is soldered to 
the circuit board at a temperature of +185°c and this tempature is 
assumed to propagate only to the lead wire during the soldering. 
Another assumption is on the deformation of the device that the 
ceramic chip and the circuit board are so much stiffer than the leads 
so all the deflections due to the thermal loading occur only to the 
leads. Finally, the device is assumed to be geometrically symmetric 
about the center so the following calculation is performed only for 
the half of the device. 

Material properties of the chip, lead and the glass-epoxy 
substrate of the circuit board are as follows; 

Thermal expansion factor of the chip, a. = 6.5 10’ 6 ppm/ c c 
Thermal expansion factor of the lead, a L = 4.47 10' 6 ppm/'c 
Thermal expansion factor of the substrate, a c = 16. 10* 6 ppm/°c 

Young's modulus of the lead wire, E = 21. 10 6 psi 

Moments of inertia of the lead wire, I = bh 3 /1 2 = 0.017 (0.005) 3 /1 2 

= 1.77 10-10 in 4 



Calculation of thermal mismatch 



At +20°c; 

Initial deformation of a lead wire after the welding 
(from +185°c to +20 = c) 

A |N|T = a AT(L) = 4.47 10 6 (-1 65^c)(0.045) = -33.2 10-6 j n 
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At +65°c; 



Deformations due to temperature changes from +20°c to 
+65°c 

A chip = a c AT(l_ c ) = 6 - 5 "'0' 6 (45 c c)(0.26/2) = 38.0 10’6 in 
A lead =ot L AT(L L ) = 4.47 1 0-6(45°c)(0.045) = 9.1 10’6 in 
A $ub =a s AT(L s ) = 16. 10‘ 6 (45°c)(0.35/2) = 126.0 10‘ 6 in 

Mismatch between the top(chip and lead) and the 
bottom portion(substrate) is 

a 65° = A top ' A bot 

( A chip + A lead + A INIT^ ’ ^ A sub^ 

= (38.0 10- 6 + 9.1 10' 6 - 33.2 10‘ 6 ) - (126.0 10’ 6 ) 
= -112.1 10-6 in 



At -20 c c ; 

Deformations due to temperature changes from +20°c to 
-20 c c 



A C hi P = a c AT ( L c) = 6 - 5 19' 6 (-40°c)(0.26/2) = -33.8 10' 6 in 

A legd = a L AT(L L ) = 4.47 10-6(-40°c)(0.045) = -8.05 10’ 6 in 

A sub =a s AT(L s ) = 16. 1 0' 6 (-40‘c)(0.35/2) = -112. 10' 6 in 

Mismatch between the top(chip and lead) and the 
bottom portion(substrate) is 

A -20" = A top ' A bot 

( A chip + A lead + A ! N IT^ ’ ( A sub^ 

= (-33.8 10-6 . 8 .1 10-6- 33.2 10' 6 ) - (-112. 10-6) 
= 36.9 10-6 in 

Since A 65 c is much greater than A. 2 q° ' n ' ts magnitude, it 
is considered that A 65 c becomes critical for the stress 
calculation. 



Analogy of temperature change effect and support settlement effect 
It can be seen that the force balance from the temperature 
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changes are the same as that from the support settlement shown in 
Figure 3. So we now treat the displacement A 65 <, from the above 

calculation as the settlement 8. 



8 = A 



65 ° 



To find the force P cause by a displacement 8 which is known 
from the previous calculation, we’ll break the frame ABC into beam 
AB and beam BC as in Figure 4 and 5. 



Denoting 



L v and 



L h as the length of the beam AB and BC, 



respectively, the tip displacements and the tip rotations at the beam 
AB due to an applied force P B and a moment Mb are given by the 

formulas 



PL. 



8 = 



3EI 



ML. 



2EI 



( 1 ) 



PL. 



0 B ~ 2EI 



ML, 

El 



( 2 ) 



and for beam BC the tip rotations can be expressed by 



ML, 
0 B = ~E\ 



(3) 



Compatibility requires that e B at the above equations (2) and 

(3) be the same so by equating these two equations an expression for 
M can be found. 



M = 



PL * 

2(L V + L h ) 

Substituting (4) into (1) 



PL. 



8 = 



L v + 4L h 



El 1 2 ( L + L h ) 



(4) 



or P can be written in terms of 8, 
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p = 



5 El 1 2(L + L h ) 



*“v + 



112.1 10-6 (21. 1Q6) (1.77 10-6) 12(0.028 + 0.045) 



(0.028) J 

Shear stress calculation 



0.028 + 4(0.045) 



= 0.0811 lb 



The shear stresses at the solder joint are obtained by dividing the force 
by footage area of the lead wire which is soldered to the circuit 
board, 



_ P _ _P_ _ 0.0811 

1 = A = bL = 0.017(0.03) 



159.0 psi 
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Analysis of unleaded device 



Figure 6 shows the geometry and dimensions of the unleaded 
device which consists of a ceramic chip fixed to a circuit board 
by a thin strip of solder. In the present application, the solder 
adhesive layer covers only a small area along the edge of the 
chip. 

Two models of the problem follows. The first model, developed 
by Suhir (1), provides an analytical determination of thermally- 
induced stresses in a tri-material assembly. The second model is 
a finite element model which utilizes a recently developed element 
which maintains inter-layer continuity along the chip-solder and 
board-solder interfaces. 

Suhir's Model 



A free body diagram of the tri-material assembly, showing the 
notation used hereafter, is shown in Figure 7. The Suhir model 
is applicable to assemblies with continuous attachment (a = I), 
as well as non-continuous attachment (a < I). The development 
is carried through by invoking equilibrium, displacement 
compatibility, and thermoelastic constitutive relations. The 

displacement compatibility relations impose the condition of inter- 
layer displacement continuity. For a detailed development of the 
analytical relations, the reader is referred to the original work. 
Here a brief summary of the development is presented, followed 
by a sequential list of the calculations that lead to the maximum 
interlayer shearing stress. 

Denoting the axial displacements at the top and bottom of the 
i th layer as u + i(x) = u'j(x) respectively, the displacement 
compatibility conditions along the interfaces of layers 1-2 and 
layers 2-3 are: 



ul ( X ) 


= u 2 (x) 


( 5 ) 


u 2 - (x) 


= u 3 * (x) 


( 6 ) 



Axial displacements u',(x), u + 2 (x), u' 2 (x), and u + 3 (x) are obtained 
by adding the individual contributions to these displacements 
arising from (a) thermal expansion, (b) axial forces, (c) shear 
forces, and (d) bending moments. The moments are related to 
the axial forces by the moment equilibrium equation which gives 
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( 7 ) 



h, + h. , h 7 +h, 

7\ (x) + -2—J. 



T 2 (x) = M 2 + M 2 + M 2 



where T, and T 2 are the axial forces on layers 1 and 2. The Tj 
axial forces are related through equilibrium to the shear stresses 
between layers, that is 



T^x) =f x x i (l)di (8) 

where ^(x) is the shear stress between the i th and (i + 1) st 
layers. The M, moments are related to the radii of curvature, p t 
through 



M^x) = - -Bl. i = 1,2,3 ( 9 ) 

P (x) 

and the Dj flexural rigidity coefficients for each layer are given 
by 



D i = ( 10 ) 

12 (1 -v^) 

where E,, h j, and are the Young’s modulus, thickness, and 
Poisson’s ratio of the i ,h layer respectively. Equation (10) assumes 
there is no variation of lateral displacement with respect to y, 
and therefore the radii of curvature of layers are all equal. 

Using equilibrium equation (7) and moment-curvature equation (9) 
yields 



l 

p (x) 



~T 2 (x) 




T 2 (x) 



\ 



D 



where 



D = D 1 + D 2 + D 1 



(ii) 



( 12 ) 



is the flexural rigidity of the assembly, and where 



^ i , i»i 






2D 



( 13 ) 



is an additional compliance due to bowing. The individual 
contributions to Eqn's (5) and (6) are the 
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• the thermal contribution is: 



a i (A r) 2 x 



(14) 



• the axial force contribution is: 

r, (£><*£ 



(15) 



where is the axial compliance of the i layer, given by 

l -v. 



k< = 

1 E.h. 



(16) 



• the shear stress contribution is: -A- | r i (x) , where k ( is the 
interfacial compliance given by 



h i 

3G, 



(17) 



• the rotational contribution is: 

if* dj 

2 Jo p(5) 



(18) 



Invoking the condition that the shear stresses be anti-symmetric 
with respect to the origin of coordinates leads to the expression 

x i (x) = C^sinh kx i = 1,2 (18) 



where k is obtained by solving a characteristic equation defined 
in Eq.(20) below. Substitution of the above expressions into 
compatibility conditions (5) and (6) yields the quartic equation 

K 12 K 33-k 4 - ( K 12*23 + *23*12) ** + 6 = 0 (20) 

where the /c^’s and A^’s are interfacial and axial compliances, 
respectively, of the tri-material assembly, given by 

K 12 = <1 + K 2 ( 21 ) 
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where 





K 23 = 


<2 + K 3 






(22) 




i° - ; 
A 23 A 2 


+ *3 + *23 






(23) 




1 0 - 1 
A 12 A i 


+ X 2 + *12 






(24) 


o 

II 


1*2 + X 2 \ 3 + X 3 X, 


+ *3*12 + *2 


(\^12 + 


V^23 ) 2 


(25) 


a) = k 2 


, Eq.(20) becomes 










Ki 2 k 33 o> 2 - (k 12 X° 3 


+ k 23 X° 2 ) o) 


+ 6=0 




(26) 


solution 


is 












<*>1.2 " 


b±y/b 2 -4ac 

2a 






(27) 


a = 


K 12 K 33 


b = - (k 12 A 


'23 + K 23 


*?2> 


(28) 


0 = 1.2) 


roots are then obtained as 










*1.2 = 


*3.4 = 


± ^2 




(29) 



The expressions presented above, valid for any tri-material 
assembly, reduce to the simpler form 

t (x) = k . Ag sinh kx (30) 

Xcoshi 



where 



k = 




K = K x + 2 1C, + K 3 



(31) 



when the Youngs modulus of the mid layer, or its thickness is 
very small, as in the case of a thin layer of solder used as an 
adhesive. In accordance with Eq.(19), the maximum shear stress 
occurs at x = f, and can be calculated from 
where 
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x = k A^Al y 

max i An 



( 32 ) 



Xmax ^ 



cosh k ( 1 -a ) 
cosh kl 



( 33 ) 



Application of Suhir Model to Unleaded Device 

Here, the closed form equations from Suhir’s model are used to 
determine the shear stresses in an uleaded device. The 
application is to unleaded device with dimensions and material 
properties as follows: 



Material #1; Ceramic chip: 



«1 = 


3.6x10® in/in/F, v 


= 0.3, E, 


= 40x10® psi, 


G 1 = 


15.4x10® psi, h, = 


.056 in. 




Material #2; 


Solder: 






«i = 


25.x10® in/in/F, v 


= 0.3, E 2 


= 1.86x10® psi 


G 2 = 


.66x10® psi, h 2 = 


.002 in. 




Material #3; 


(board): 






« 3 = 


16.x10® in/in/F, v 


= 0.4 E 3 


= 2.0x10® psi 


g 3 = 


.77x10® psi, h 3 = 


.056 in. 


* 


above parameters lead to the following compliances: 


Chip: 


D, = 643 lb. in., 


. ^ = 


.3125x10® in. /lb 




tc, = 1.212x10'® 


in. /lb. 




Solder 


: D 2 = .0231 lb. in, 


., A 2 = 64.52x10® in. /lb., 




k 2 = 2.525x10® 


in. /lb. 




Board: 


D 3 = 32.2 lb. in. 


. a 3 = 


6.25x10® in. /lb., 




/c 3 = 24.24x10® 


in. /lb. 





The other parameter values necessary for a solution are: 

A 12 = 1 .378x1 0 6 in. /lb. A 23 = 1.378x1 O' 6 in. /lb. 
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3.737x10® in. /lb. 



26.77x10® in. /lb. 



A 12 ’ = 64.83x1 O' 6 in. /lb. 



A 23 ' = 70.77x10® in. /lb. 



A 12 ° = 66.21x10'® in. /lb. 



Xj = 72.15x10® in. /lb. 



Substitution of these parameters into Eq.(25) yields 

b = .7897x10® (in. /lb.) 2 



which upon substitution into Eq.(26) yields k = 19.85 in.' 1 . 
Finally, the maximum shear stress at the solder interface is 
obtained from Eq.(32) as 



The simpler approximate expression for k from Eq.(31) yields k 

= 19.61 in.' 1 , which, in turn, gives r max = 4310 psi. 

Finite Element Model 

An alternate model of the tri-assembly problem using a finite 
element method was undertaken. Continuity along the interfaces 
of adjacent layers is maintained by selection of the six degree 
of freedom element shown in Figure 8. The six degrees of 

freedom are as follows: 

b , is the axial displacement at the lower left corner 

b 2 is the axial displacement at the upper right corner 

b 3 is the lateral displacement at the left end 

b A is the axial displacement at the lower right corner 

b 5 is the axial displacement at the upper right corner 

b 6 is the lateral displacement at the right end 

The axial displacement field u(x,y) is taken as linear interpolation 
in both the axial, x, and transverse, y, directions. That is, 



r, 



max 



4360 psi 
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and the lateral (i.e. .transverse) displacement field is 
function of x, that is, 

v(x) = N 2 (x) 63 + n 2 ( x ) 5 6 

where the linear interpolation functions are: 

H i (y) = \ ^ h ~y ) 






y 

h 



N x (x) = 1 (1-x) 



N 2 (x) = ^ 



The strain-displacement relations are: 






dv A7, (•*) 



dx 



/s s x (x) 

( 5 2 _ 5 i ) + — S; — <» 5 -» 4 ) + 



On the element level, we have the bending effects 



Mb = f b 



where 



and 



6 b — < 5 ^ 64 ^5^ 



K b = f f BEB T dydx 
J 0 Jo 



a linear 

( 35 ) 

( 36 ) 

( 37 ) 

( 38 ) 

( 39 ) 

( 40 ) 
^ ( 41 ) 

( 42 ) 

( 43 ) 

( 44 ) 
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and 



and 



r dN. 

BT - ( ^ 



dN, 

id H > 



dN, 

-d H ' 



^2 

dx 



H,> 



Fb = 1 1 1 h BEB T dydx 



which gives 



'2 l-2-i 

Eh 1 2 - 1-2 

B 6 1 - 2-121 
1-212 



Similarly for the shear effects 

Ms = o 



where 

= < 6 X b 2 6 3 6 4 6 S S 6 > 

and 

K s = f 2 f h B'GB ,T dydx 
J 0 J 0 

where 



, T dH , 

B /r = <w.M 

ay 






a^ 

ay 



awi w <^1 

"a37 2 “ay 





which gives 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 



(51) 
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-a 



-a 



G _ G 
2 2 



a -a 



a -a 



a 



a 




( 52 ) 





The bending and shear element stiffness matrices are combined 
to form the element stiffness matrix, which in turn are again 
combined to form the system stiffness matrix. In accordance with 
standard FEM techniques, the system matrix equations are solved 
for the system displacement degrees of freedom. The element 
strains are then calculated from Eq.(40) and (41) above. Finally, 
the average bending stresses at the bottom and top of an 
element are obtained from the stress-strain relations 



Similarly, the average element shear stresses are calculated from 
the strress-strain relations 




( 53 ) 



and 




( 54 ) 




( 55 ) 



and 




( 56 ) 
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The finite element mesh, shown in Figure 9, consisting of 104 
elements and 149 degrees of freedom produced an interface 
maximum shear stress of 3650 psi. This FEM calculated stress is 
at the midpoint of the right most element. This stress, obtained 
from an FEM analysis is at the midpoint of the interface at the 
right most element. Quadratic extrapolation to approximate the 
hyperbolic sine expression for the r shear stress brings the FEM 
predicted value to about 3800 psi. Thus the predicted values of 
the maximum shear stress from the Suhir model and the FEM 
model differ by less than 15 percent. 
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Conclusions 



This investigation has shown that under the same thermal 
environment, the stiff unleaded solder joint has significantly higher 
thermoelastic stresses than the flexible leaded solder joint. Since 
the unleaded joint has a stronger suppression of displacement 

than the leaded joint, and thermal stresses result from 

displacement constraint, the result is not an unexpected one. In 
the particular application presented in the previous sections, the 
stresses in the unleaded connection were twenty times the 
stresses in the leaded connection. This large degree of difference 
between the stresses in the two cases is surprising. 

As previously noted, the comparison is for the case where both 
connections are subjected to the same thermal environment. This 
assumption does not take into account the fact that the thermal 
field is very likely to be different for the two types of joints. 

One might expect that the leaded joint with the ’separated’ 

components may have larger thermal gradients than the unleaded 
joint where the components are in close proximity to each other 
and therefore may have a more uniform thermal field. In fact, the 
solution to the stress problem really requires its coupling to the 
thermal (energy balance) problem. This later coupled (and 

transient) problem is presently under investigation as a NPS 
thesis topic under the direction of D. Salinas and Y.W. Kwon (2). 
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Figure 1. Geometry of leaded device 
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Figure 2. Cross section of lead wire 
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Figure 3. Support movement 6 
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Figure 4. Deformation of beam AB 
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Figure 5. Deformation of beam BC 
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Figure 6. Unleaded device 
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Figure 7. Force and moment equilibrium of components 
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Figure 8. Six degree of freedom element 




Figure 9. 



Finite element mesh 
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